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POLYPHASE FILTER WITH STACK SHIFT CAPABILITY 

BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

The present invention relates to a discrete Fourier transform (DFT) filter bank, 
and more particularly to a polyphase filter bank v^th stack shift capability which can be 
implemented with significantly less hardware than known filter banks. 



2. Description of the Prior Art 

Digital filter banks are generally known in the art. Examples of such digital 
filter banks are disclosed in U.S. Patent Nos. 4,107,470; 4,393,456; 4,792,943; 4,839,889; 
5,436,940; and 5,606,575, hereby incorporated by reference. Such digital filter banks are 
implemented with digital signal processors and are used in a wide variety of applications 
including telecommunications applications, for example, cellular telephone applications. 

Such digital filter banks are known to be used as spectrum analyzers as 
illustrated in FIG. lA and fi-equency synthesizers, as illustrated in FIG. IB as discussed in 
detail in Multi-Rate Digital Signal Processing , by Ronald E. Crochiere and Lawrence R. 
Rabiner, Prentice Hall, copyright 1983, hereby incorporated by reference. As is generally 
known in the art, the principles relating to frequency synthesizers are normally applicable to 
spectrum analyzers and vice versa. 

Referring to FIG. 1 A, a broad-band input signal x(n) is applied to a digital 
filter bank, configured as a spectrum analyzer, which divides the input signal x(n) into K 
frequency channels, 0 through K-1, which correspond to K output signals x^Cm), 
x,(m)...Xk(m)...XK.,(m). 
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Similarly, the system illustrated in FIG. IB may be used to synthesize a 
plurality of input signals Xo(m), Xi(m)...Xk(m)... XK.i(m), defining K input channels, 0 through 
K-1. In this embodiment, K input signals Xo(m), x,(m)...XK.i(m) representing K channels 0 
through K-1 are synthesized by the digital filter bank to form a broad-band output signal 
5 x(n). 

FIGs. 2A and 2B represent models for a DFT filter bank analyzer and a DFT 
filter bank synthesizer, respectively, as set forth in Multi-Rate Digital Signal Processing . 
supra. For simplicity, only a single channel is illustrated. As illustrated in FIGs. 2A and 
23, the sampling fi-equency of the analyzer is decimated or decreased while the sampling 
1 0 firequency for the synthesizer is interpolated or increased. As such, such DFT digital filter 
banks are multi-rate in nature. 

As shown in FIGs. 2A and 2B, the DFT filter banks are modeled by a low pass 
filter h(n), f(n) and a complex modulator. In addition, as discussed above, the sampling 
frequency is either decimated by a factor M as in the case of a filter bank analyzer or 
p 15 interpolated by a factor M in the case of a filter bank synthesizer. As shown in FIG. 2 A, in 
^ the case of the analyzer, an input signal x(n) is modulated by a complex modulator of the 

J: form e"^"'*" and low pass filtered by a filter h(n). The sampling rate is then reduced or 

' decimated by a factor M to generate the channel signals Xk(m). The filter h(n) determines 

^ the width and the frequency response of each of the channels. The filters h(n) for each 

i 20 channel are identical. 

0 For a DFT filter bank with uniformly spaced filters and even type stacking 

arrangements, the center fi-equency of each of the channels is given by equation 1 below: 



27tk 

m = ^.k = 0.1....,K-l 



(1) 



By defining WK=e ^^^"^^ the complex modulation function can be rewritten as shown in 
2 5 equation 2 below: 



(2) 
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The channel signals may then be expressed as set forth in equation 3: 



00 



XK(m)=^h(mM-n)x(n)wt,lc = 0,l,...,K-l 



(3) 



With respect to the DFT filter bank synthesizer, for example as shown in FIG. 



2B, all of the input channel signals Xo(m), x,(m)...Xk(m)...XK.,(m) are interpolated to a higher 
sampling rate and modulated back to the original spectral location. The synthesizer then 
adds all of the channel signals together to produce a single output signal x(n). Each of the 
input channel signals Xo(m), x,(m)...x,,(m)...XK.i(m) are interpolated by a factor M and filtered 
with an interpolation filter f(n) and modulated by a complex modulation fimction Wj^^"*^ = 
^-Jmkn gjjjf^ ^YiQ charmel signal back to its original location co^. Equation 4 represents the 
output of each channel signal: 



Since the synthesizer output signal is the sum of all of the channel signals as shown in 
Equation 5, Equation 4 can be rewritten as shown in Equation 6: 



00 



x,(n) = w'^ 2 X,(m)f(n-mM),k = 0,l,...,K.l 



(4) 



x(n) = ^f^h(n) 



(5) 



00 1 K-I 



(6) 
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In order to improve the efficiency of such DFT filter banks, such DFT filter 
banks are known to be implemented as polyphase filters. More particularly, the decimator 
M and low pass filter h(n) implementation of the fi:equency analyzer and the interpolator M 
and filter f(n) of the frequency synthesizer are replaced with polyphase filter structures Pp 
5 (m) as illustrated in FIGs. 3 A and 3B, respectively. As mentioned above, both the analyzer 
and the synthesizer are multirate devices. The complex modulation as well as the polyphase 
filters Pp essentially shift the frequency of the prototype filters (h(n) or f(n)), for example as 
shown in FIG. 4. As shown, the center frequency of each of the M frequency channels is a 
multiple of the sampling frequency Fg, divided by the number of channels. As shown in 
1 0 FIG. 4, the transition bands of the low pass filter get shifted and become gaps in the 

frequency response of a DFT filter bank. In order to resolve this problem of gaps in the 
frequency response, DFT filter banks are utilized with stack shift capability. An example of 
such a DFT filter bank is illustrated in FIG. 5 A, a polyphase filter bank with M channels, 
each channel including the following twiddle factors: tap twiddle, DFT twiddle and output 
1 5 twiddle, as discussed in detail Multi-Rate Digital Signal Processing , supra. Each tap 

twiddle factor includes a complex term in the form^^^''*' , wherein k^ is the stack shift 
\f[ factor. The tap twiddle is used to provide the stack shift. Unfortunately, such a DFT filter 

bank with tap twiddle requires an inordinate amount of hardware to implement. For 

example, to implement an 8 phase 203 tap filter as illustrated in FIG. 5 A, 21 1 selectable 

1 1 

£ 20 negators are required to produce even/odd stack shifts. Unfortunately, this amount of 
:^ hardware duly complicates the system and raises the cost. Thus, there is a need for a more 

hardware efficient implementation of a digital filter bank with stack shift capability. 

SUMMARY OF THE INVENTION 
2 5 Briefly, the present invention relates to a discrete Fourier transform (DFT) 

filter bank with stack shift capability, configured, for example, as an analyzer, but equally 
applicable to a synthesizer, and implemented as a polyphase filter. The DFT filter bank 
includes a complex modulator for multiplying the input signals x^{rci)^ Xi(m)...x^,.,(m) by a 

complex modulator ^ (m). In order to keep the coefficients for the polyphase filter real, 

30 the coefficients are modulated by a DFT twiddle factor e^^""^"^ , xhe use of the DFT 
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twiddle factor enables the hardware to be reduced significantly. For example, for an 8 phase 
203 tap filter, only 8 selectable negators are required as opposed to 21 1 for known DFT filter 
banks. 




5 DESCRIPTION OF THE DRAWINGS 

These and other advantages of the present invention will be readily understood 
with reference to the following specification and attached drawing wherein: 

FIG. 1 A is a generalized block diagram for a K channel filter bank analyzer. 
FIG. IB is similar to Fig. la but for a K channel filter bank synthesizer. 
1 0 FIGs. 2A and 2B are block diagrams showing additional details of the DFT 

filter bank analyzer and synthesizer illustrated in FIGs. la and lb. 

FIGs. 3A and 33 illustrate commutator models of a DFT filter bank analyzer 
and synthesizer implemented as polyphase filter networks. 

FIG. 4 is a graphical illustration of the frequency response of a known DFT 
1 5 filter bank analyzer or synthesizer. 

FIG. 5A is a block diagram of a conventional DFT polyphase filter bank 

analyzer. 

FIG. 5B is a block diagram of a DFT polyphase filter bank analyzer in 
accordance with the present invention. 
2 0 FIG. 6 is a block diagram of an input twiddle block in accordance with the 

present invention. 

FIG. 7A is a block diagram of a polyphase filter structure for use in accordance 
with the present invention. 

FIG. 7B is a block diagram illustrating the use of multiple compensation 
2 5 vectors for use with the polyphase filter illustrated in FIG. 7A. 



DETAILED DESCRIPTION 
The present invention relates to a polyphase DFT filter bank with stack shift 
capability. Although the principles of the invention are described and illustrated herein with 
3 0 respect to a frequency analyzer, as is generally known in the art, such principles are also 
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applicable to DFT filter banks implemented as fi'equency synthesizers. All such 
embodiments are intended to be covered by the scope of the appended claims. 

As mentioned above, the DFT filter in accordance with the present invention is 
adapted to significantly reduce the hardware for implementing DFT filter banks 
5 implemented as polyphase filter networks with even/odd stack shift capability. As discussed 
above, for an 8 phase 203 tap filter, known filter banks require 211 selectable negators while 
only 8 are required for the present invention. 

As discussed in Multi-Rate Digital Signal Processing , supra, a stack shift may 
be implemented by shifting the filter response in mathematical terms. More particularly, as 
1 0 illustrated in FIGs. 5 A and 5B, the DFT filter bank analyzer includes M polyphase filters 

Po(m)...pp(m).,.pj^.i(m). In order to accomplish the odd/even stack shift, the filter response of 
each of the polyphase filters po(m)...pp(m)...pj^.i(m) is multiplied by a term Wm^ = ^j^^^o^^ 
resulting in all of the polyphase filters Po(m)...Pp(m)...P^^_i(m) becoming generally complex 
filters. In addition, the output of each of the polyphase filters is modulated by a factor 

15 ^^^'^^"^^^ , where ko = 0 or 0.5. The output sum is multiplied by or (-1)"". With such 
a configuration, the foUov^ng additional hardware is required to implement the even/odd 
stack shift: selectable negation for all filter taps; selectable negation for the output signal; 
and selectable DFT twiddle factors. 

2 0 With such an implementation, the output signal is given by equation 7 below: 

Xp(m) = x(mM^ p),p^O,l...,M ' 
Pp(m)^h(mM-p),p = OJ M- 



p-O 



(7) 



2 5 where m is the time index, p is the input/polyphase filter branch (from 0 to M-1), k is the kth 
channel of channelizer (fi-om 0 to M-1), ko denotes the channelizer stacking mechanism (0 
for even stacking, 0.5 for odd stacking), and M is the total number of branches/channels. 

Xp(m) are the M input branches, ^ is the filter tap modifying factor, Pp(m) are the 8 
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polyphase filter branches, ^ 



j2n/MkoP 



M-l 



is the DFT twiddle factor, and S ^'^"^^ denotes the 

p=0 



DFT operation. 

In order to reduce the amount of hardware needed to implement the filter taps, 
the DFT filter bank in accordance with the presenting invention eliminates the need to 
5 selectively negate all filter taps and the output signal and instead selectively negates the 
input signal only, thus significantly reducing the amount of hardware required to implement 
the filter bank. Referring to FIG. 5B, the DFT filter bank analyzer in accordance with the 
present invention is illustrated and identified with the reference numeral 20. As shown, the 
DFT filter bank 20 is adapted to analyze input signals Xo(m)...Xp(m)...Xi^.,(m), which are 
1 0 polyphase components of X(n), and provide output signal. In order to shift the input 

fi-equencies by co, the input signals Xo(m),..Xp(m)...XM.i(m) are multiplied by a factor e"^"^", 
where n is the time index of the input signal. This operation results in a generally complex 
signal, which, in turn, requires a complex filter to process and thus doubles the size of the 
filter. 

1 5 The filters are implemented as polyphase filters 

Po(m)...pp(m)...pM_i(m). Rather than modify the frequency response of each of the filters in 
the polyphase filter bank as is known, the present invention shifts the frequency of the M 
input signals Xo(m)...x„(m)...XM.i(m), by co where co=7r/M and M equals the total number of 
channels. The multiplication factors preceding each branch or channel of the polyphase 

2 0 filter bank take on the following form: 



25 



Branch No. 
0 
1 

M-2 
M-l 



First M Input Samples Second M Input Samples 

^(ji/M)*0 ^(jtM)*M 

^(7t/M)* 1 ^(jt/M)*(M+l )_^^(7i/M)* 1 



^(7i/M)*(M-2) 
^(Ji/M)*(M-1) 



^(7c/M)*(2M-2)__^(7:/M)*(M-2) 
^(;iflM)*(2M-l)__^(7i/M)*(M-I) 



Third M Input Samples 

^(7i/M)«(2M)^^(K/M)»0 
^(7u/M)«(2M+l)_^(7:/M)*l 

^(ji/M)*(3M-2)_^(7i/M)*(M-2) 
^(7:/M)*(3M-l)_^(jc/M)*(M-l) 



The fourth M input samples have multiplicative factors identical to those of the 
30 second M input samples and so on. As can be observed above, the multiplicative factors for a 
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particular branch have a repetitive pattern where every other one has the same value and 
adjacent factors are negatives of each other. As such, these multiplicative factors may be 
implemented with relatively simple hardware by the input twiddle and DFT twiddle circuit 
illustrated in FIG. 5B. 

5 FIG. 6 is a block diagram of the input twiddle circuit. The input twiddle circuit 

provides a ±1 factor for the multiplicative factors as discussed above. In particular, the input 
twiddle circuit may be implemented as a selectable alternating inverter. The input twiddle 
circuit includes a divide by two counter which essentially divides a clock input by two. The 
output of the divide by two counter 22 is ANDed with a stack shift factor ko by way of an AND 

1 0 gate 24, used to generate the least significant bit output Isbout. The output of the AND gate 24 
also drives a multiplexer 26 in which the input signal XQ(m)... Xp(m)...XM.i(m)... and its 
complement are applied by way of an inverter 28 to multiplexer 26. The output of the 
multiplexer 26 is an output signal OUT that is applied to the input polyphase filters Po...Pm.i. 
In FIG 6, kfl is set to logic 0 when the filter is to be configured as an even-stack 

15 filter (ko =0 in equations above). It is set to logic 1 when the filter is to be configured as an 
odd-stack filter (k^ = 0.5 in equations above). Thus when the filter is configured as an even- 
stack filter, the ouput of the AND gate 24 is always logic 0 in which case the MUX 26 selects 
the IN signal. In other words, the IN signal passes through the circuit imaltered. When the 
filter is configured as an odd-stack filter, the output of the AND gate 24 is the same as the 

2 0 ouput of the clock divide block 22 (since k^ is logic 1), which is alternating high and low. In 
this configuration, the MUX 26 ouputs the signals IN and IN inverted (output of the INV 28) 
on alternating clock cycles, effectively negating the input on every other clock cycle. 

Negation in two*s complement of a number is accomplished by inverting all bits and 
adding 1 to the LSB. Referring to FIG 6, if the Isbout and OUT signals are added together, 

2 5 then input negation on alternating clock cycles is complete. However, the addition of the 

Isbout and OUT signals needs to be propagated all the way up to the MSB of the sum; 
potentially a slow operation. In order to avoid this operation, compensation vectors are used. 
Such compensation vectors are known in prior art, for example as disclosed in, "A Silicon 
Complier for High-Speed CMOS Multirate FIR Digital Filters," Robert Hawley M.S. thesis, 

3 0 University of California, Los Angeles, 1 99 1 , hereby incorporated by reference. With the LSB 

addition incorporated into the compensation vectors, only an inverter is needed to implement 
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the input twiddle circuit. This eliminates the need to do a true two*s complement negation and 
the carry propagation associated with it, allowing high speed operation of the filter without 
additional pipelining. 

As mentioned above, the input twiddle circuit provides an alternating 0 and 1 output 
5 at the Isbout output. This least significant bit Isbout signal is applied to the polyphase filters Pp 
to complete the two's complement negation through the use of compensation factors as 
discussed below. 

FIG 7 illustrates one branch of a polyphase filter implemented in CSD coefficients, 
for example, as disclosed in the Hawley reference, supra. Normally, only one compensation 

1 0 vector is required per filter branch. However, in order to implement a stack shift, one 

compensation vector is required for an even stack (ko=0) and two additional compensation 
vectors are required for an odd stack (ko=0.5). The compensation vector for the even stack 
configuration can be derived using the method described in Hawley, supra. The compensation 
vectors for the odd stack configuration can be derived similarly, with modified terms for LSB 

1 5 addition due to 2's complement negation of CSD coefficients. 

As disclosed in the Hawley reference, supra, an LSB addition term is added to the 
compensation vector when a CSD coefficient of a filter tap is negative. To illustrate, let C be a 
positive number and X be the input signal. An LSB term is added to the compensation vector 
if a CSD coefficient is negative (-C*X), No additional term is required in the compensation 

2 0 vector when a CSD coefficient is positive (C*X). But since -C*(-X)=C*X and C*(-X)=-C*X, 
when the input signal is negated (-X) by the input twiddle, an LSB term is required in the 
compensation vector if a CSD coefficient is positive (C*(-X)=-C*X), while no additional term 
is required if a CSD coefficient is negative (-C*(-X)=C*X). 

For an odd stack shift (ko=0.5) the inputs Xo, Xj, Xj, X3.. .become x^,, -Xi, Xj, -X3..., the 

2 5 filter outputs for a single phase in terms of the taps ho, hj, hj, ha.. .become yo = h^, * Xq, yj = ho * 
(-xO + h, * (xo); y2 = ho * (x^) + h,* (-x,) + h^ * (x^); y3 = ho * (-X3) + h^ (x^) + h^* (-x,) + ha * Xq. 
These filter outputs are completed by two compensation vectors identified in Fig 7, for 
example, as comp vector 1, comp vector 2. The first compensation vector comp 1 vector is 
used to add an LSB to negative CSD bits of even taps (ho, h2, h4...) and positive CSD bits and 

30 odd taps (hj, h^, hy.,) and is applicable for outputs yo, y2... y2n- The second compensation vector 
comp 2 vector is used to add LSB to negative CSD bits of odd taps and positive CSD bits of 
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even taps for outputs Yi, ya— y2n+i- 

These different comp vectors, comp vector 1, comp vector 2, are applied to a MUX 
32 which, in turn, is under the control of the signal Isbout (FIG 6). The MUX 32 is applied to 
one input of a MUX 34. The MUX 34 is used for even/odd stack shift. In particular, for even 
stack shift, an even stack compensation vector is applied to the polyphase filter structure. 
During conditions of an odd stack shift, one of the odd stack shift compensation vectors, comp 
vector 1, or comp vector 2 is applied to the MUX 34 under the control of the signal Isbout. 
As mentioned above for an odd stack shift, the signal Isbout alternates between 0 and 1 logic 
states, thus providing the proper compensation factor for the particular sample period. 

In order to maintain the coefficients for the polyphase prototype filter real, the 
complex portion of the multiplicative factors is implemented at the end of each branch, shown 

as DFT twiddle factors . The filter bank 20 in accordance with the present invention 

using the stack shift scheme discussed above is thus described by equation (8) below: 

Y, (m) = |]e>2«*p/A/ ^ ^-Jl^oPfM ^ ^-72;*o. * ^^^1 ^ ^ 0,1,..., M - 1 

p=0 

(m) = x(Mm -p\p = 0,1,...,M - 1 
T^(m) = h(Mm + pXp = 0,1,...,M - 1 

(8) 

where m is the time index, p is the input/polyphase filter branch (fi-om 0 to M-1), k is the kth 
charmel of channelizer (from 0 to M-1), ko denotes the channelizer stacking mechanism (0 for 
even stacking, 0.5 for odd stacking), and M is the total number of branches/channels. Xp(m) 
are the M input branches, e'-^^""^"'" is the input twiddle factor, (m) are the 8 polyphase filter 

branches e'-'^'^^''^ is the DFT twiddle factor, the DFT operation is denoted by equation (9). 

M-I 

(9) 
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Obviously, many modifications and variations of the present invention are possible 
in light of the above teachings. Thus, it is to be understood that, within the scope of the 
appended claims, the invention may be practiced otherwise than as specifically described 
above. 

What is claimed and desired to be covered by a Letters Patent is as follows: 




